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Abstract 

The recently observed quantum states of neutrons bound in a gravitational field are analyzed 
in the framework of one-parameter isospectral hamiltonians. Potentials isospectral to the usual 
Newton potential are explicitly constructed for the first time, then constrained using measured 
properties of the neutron gravitational states. The corresponding wave functions and the neutron 
fluxes are also calculated and analyzed in a simple model, including the ground state and the 
excited state contributions. The constructed isospectral potentials are discussed as candidates 
for a possible modification of Newton's law at a submillimetre scale. Our results indicate that 
significant deviations from the Newtonian gravity at submillimetre distances could be compatible 
with experiment. 
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I. INTRODUCTION 



In recent years there has been a significant progress in the study of gravity on a submil- 
limetre scale. It has been precipitated by the first measurements [3, 0] of the quantum levels 
of ultra-cold neutrons propagating under gravity within a submillimetre-wide slit between a 
bottom mirror and a top absorber. These results have profound implications as they could 
put experimental constraints on proposed extensions of the Standard Model (see e. g. jl, 0, 0] 
and references cited therein). 

The data on the quantum gravitational states of neutrons have been analyzed in several 
phenomenological models [l|, H, @, 0, @, @, E3] with the aim of exploring possible deviations 
from Newton's law at a submillimetre scale. Although the data can be explained in models 
with unmodified gravity and a more-or-less sophisticated description of the neutron guide 
system 0, @, 0, [l0| , the current experimental uncertainties do leave some room for possible 
modifications of Newton's law at submillimetre distances jl, @]. 

In this paper we consider this problem from the point of view of quantum-mechanical 
isospectral hamiltonians. Construction methods and various applications of isospectral 



hamiltonians can be found in Refs. [11|, [12|, [13j and in monographs [14], 1151 ] . An important 



feature of isospectral hamiltonians is that despite having quite distinct potentials (which we 
will call "isospectral potentials"), they all have identical energy spectra. This remarkable 
property lies at the core of our motivation to explore possible deviations from the Newton law 
by studying potentials isospectral to the usual gravitational potential. In this way we ensure 
that any possible modification of the Newton potential preserves the quantum-mechanical 
energy spectrum and is thus compatible with the measured properties of the gravitational 
levels of neutrons. 

We will first construct a one-parameter family of quantum-mechanical hamiltonians 
isospectral to the hamiltonian with the usual Newtonian potential of a neutron in the grav- 
itational field, and calculate the corresponding wave functions of the ground and excited 
states. We will use measured properties of the neutron gravitational quantum states [H, 0] 
as experimental constraints to select only those isospectral hamiltonians which are compat- 
ible with these experimental data. Specifically, the characteristic step-like dependence of 
the measured neutron flux on the width of the slit will be used to select the most likely 
isospectral potentials. We will explore the possibility that these isospectral potentials could 
describe deviations from the Newtonian gravitational law at a submillimetre scale. (There 
are important qualifications which limit the conclusions drawn from a comparison of the 
present simplified model with the existing data; this will be explained in the text, and we 
will suggest a modified experimental arrangement to more directly test our results.) The 
results of our model analysis indicate that the experimental data allow for appreciable mod- 
ifications of Newton's gravitational potential at heights of the order of tens of micrometres. 

II. ONE-PARAMETER ISOSPECTRAL GRAVITATIONAL HAMILTONIANS 

We consider a neutron of mass m, at a height z over the surface of a perfect mirror. Such 
a system is described by a one-dimensional hamiltonian 
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with the potential 

V(z > 0) = mgz, V{z < 0) = oo , (2) 

where g is the usual gravitational acceleration. 

The solution of this quantum-mechanical problem is well known [lqj . The neutron en- 
ergy levels E n are proportional to the zeros of the Airy function Ai [171 ]. Ai ~ 2.34, A2 ~ 
4.09, A 3 « 5.52 ... (Ai(-X n ) = 0): 

E n = \ n [mg 2 h 2 /2] 1/3 , n = 1, 2, 3, . . . , (3) 
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or numerically: Ex « 1.41 peV « 2.25 x 10~ 24 erg, £ 2 « 2.46 peV « 3.92 x l(T 24 erg, £7; 
3.32 peV ~ 5.29 x 10~ 24 erg, . . . These energy levels can be also written with high precision 
in the WKB approximation; see [l6| for details. The neutron wave functions are 



> 0) = C n Ai (~ - A n ) , < 0) = , (4) 

where I = [h 2 /(2m 2 g)} 1 ^ 3 ~ 5.87/im is the characteristic length and the constants C n are 
determined from the normalization condition dzijj 2 (z) = 1. 

We are interested in a family of hamiltonians, parametrized by a real number p, 

which have the same energy spectrum E n as the original hamiltonian H(z). Many deriva- 



tions of such isospectral hamiltonians have been developed in the past (see [111 . Il2l . 1 1 31 ] and 



references therein). One modern approach is usually formulated in the language of super- 



symmetric quantum mechanics [13] , 1 151 ] . Since the focus of this paper is the problem of 
neutron gravitational states, rather than isospectral hamiltonians in general, we will simply 
use the necessary results from [15] without proof. 

Given the original potential V(z) and the normalized ground state wave function ipi{ z ) = 
ip(z), the isospectral potential in Eq. (j5J) can be written as 

V(z;p) = V(z) - ^--^[\n(l(z)+p)} , (6) 

where 

z 

X{z)= [dz>^ 2 (z>). (7) 



Since ip{z) is normalized to unity, it follows from Eq. (JTj) that T(z) G [0, 1] and therefore, to 
ensure that V(z;p) in Eq. (jBJ) is finite, p ^ [—1,0]. 

Several members of the one-parameter family of potentials V(z; p) isospectral to the linear 
Newtonian potential mgz are shown in Fig.[H demonstrating the deformation of the potential 
with varying parameter p. (The figure will be further explained and discussed in the next 
section.) To our knowledge, this is the first time that these isospectral potentials have been 
explicitly calculated. Being isospectral to the exactly solvable quantum-mechanical system 
with the potential mgz, each of the potentials in Fig. [1] also furnishes an exactly solvable 
system [l3(, which makes them an interesting object for future studies in their own right. 
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FIG. 1: Examples of isospectral potentials of a neutron, Eq. ([6]), for p > (left panel) and p < — 1 
(right panel). The energies of the ground and first excited states E\ = 2.25 x 10~ 24 erg and 
Ei = 3.93 x 10 -24 erg are shown by the horizontal lines, and the diagonal dotted line is the 
Newtonian potential mgz. The experimental uncertainties of the measured ground state turning 
point z e ^ Lp G [10.3, 14.1] /im and of the excited state z%° p G [19.2, 24.0] fim are indicated by the 
intervals between the squares on the corresponding energy levels. The potentials compatible with 
the measured turning point cross the horizontal line within these interval. 



The ground state wave function corresponding to the isospectral potential Eq. <^ can be 
written 



i/)i(z;p) = C(p) : 



(8) 



I(z)+p' 

as can be checked by substitution into the Schrodinger equation. The proportionality con 
stant C(p) can be obtained from the normalization condition J^°dzipf(z;p) = 1: 



\C(p)\ 



jl? z 


\ M*) 1 


I(z) +p_ 





(9) 



Note that the wave function changes with changing p, while nevertheless corresponding to 
the same ground state energy level E\. The original potential and ground state wave function 
are restored from their isospectral counterparts in the limits p —>■ ±oo. 

All isospectral potentials V(z; p) have the same full spectrum of energies as the original 
Newtonian potential V(z), including the ground state and the excited states. By con- 
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trast, Eq. (jSj) is valid for the ground state wave function ipi(z;p) = ip(z;p) only. The 
normalized wave functions ip n +i(z;p) of the excited states of the isospectral hamiltonians 
can also be constructed [l3| from their unmodified counterparts ip n+ i{z) (n 



1,2,3,...): 



ip n+1 (z;p) = V>n+l(z) + 



h 2 di{z) 

2mE n+1 [T(z) + p] dz 



d 

dz 



1 dipi(z) 
ipi(z) dz 



(io) 



where the energies E n+1 and the integral 1{z) are given in Eqs. ^ and ([7j), respectively. 



III. EXPERIMENTAL CONSTRAINS ON THE ISOSPECTRAL POTENTIALS 

In this section we will determine the values of p for which the isospectral potentials 
V(z\ p) are compatible with the experimentally observed [l|, 0] neutron gravitational states. 
The lowest height at which the flux of ultra-cold neutrons has a "step", corresponding to 
the lowest quantum state, has been measured to be z\ v = (12.2 ± 1.8 syst ± 0.7 siai )/zm, 
with the estimated total error Azl xp « 1.9 fim. The (much less pronounced) next leveling 
of the flux, corresponding to the first excited state, has been extracted to occur at z%" p = 
(21.2 ± 2.2 syst ± 0.7 s t a t)/im, yielding a total error Az^ v ~ 2.3 fim. Within the error bars, 
these heights are consistent with the theoretical turning points of the ground and first 
excited states associated with the Newton potential mgz (e. g., see Refs. jg, Gil). It should 
be emphasized, however, that Z2 has so far been measured with much worse resolution than 
Z\ 0]; therefore, in our model, we will treat the constraints stemming from z^ v as less firm 
than those from z exv . 

Our objective is to select those isospectral potentials V(z;p) for which the semiclas- 
sical turning points of a neutron with the energy E± t 2, obtained from the equations 
V{z\ t 2\p) = Ei t 2, lie within the experimental error bars, i. e. z\ v G [10.3, 14.1] fim and 
z exp ^ [10.2, 24.0] fim. It can be seen from Fig. [T]that the unmodified Newtonian potential 
mgz itself satisfies this condition, as it should. Among the isospectral potentials, only those 
with p > 1.08 and p G [—1.13, —1.04] are compatible with the measured ground-state turn- 
ing point z\ v , i. e. they cross the ground state energy level between the squares indicating 
the experimental uncertainty interval. The error bar associated with z% xp does not put addi- 
tional constraints on the potentials with p > 0, but for p < —1 it appears to marginally allow 
only the potential with p = —1.13. However, we should reiterate here that the constraints 
stemming from z e 2 xv are at present much less reliable than those from z\ v . (The potentials 
with p < —20 also fulfill our selection criteria, but their deviations from the unmodified mgz 
potential are extremely small; therefore, we do not show them and will not discuss them 
further.) 

It is interesting to note that the ground state wave functions corresponding to the 
range p > 1.08 have a behaviour qualitatively different from those corresponding to 
p G [—1.13, —1.04]. In particular, the latter wave functions describe neutrons predominantly 
localized in a potential well above the turning point, which is quite distinct from the wave 
function for the usual Newton potential, or the isospectral potentials with p > 0. In the next 
section we will exploit the properties of the ground state wave functions to put additional 
phenomenological constraints on the potentials with p > 1.08 and p G [—1.13, —1.04], and 
in Section |V] we will consider the contributions from the first excited state for the potentials 
with p = 1.08 and p = —1.13. 
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IV. GROUND STATE NEUTRON FLUX 

The experimental method [l| is based on a measurement of the step-like shape of the flux 
of the neutrons passing through a slit between the bottom mirror and the top absorber. The 
measured neutron flux is typically analyzed in detailed models including the interaction of 
the neutrons with the absorber, which is related to finite lifetimes of the quantum states. 
Such analyzes were developed in Refs. 0, 0], where no modifications to Newton's potential 
were considered, while some assumptions were made about the widths of the quantum states, 
and the neutron energies and populations of the states were fitted to reproduce the measured 
neutron flux as a function of the slit height. A careful modeling of the absorber is crucial in 
these calculations to achieve a good description of the data. 

By contrast, the focus of the present model is to explore the possibility that isospectral 
potentials could describe small-range deviations from the Newton gravitational potential. In 
this section we study whether it might be possible to use experimental data to put additional 
constraints on the allowed isospectral potentials obtained above, all of which are compatible 
with the experimentally observed ground state as described above. (The contributions of the 
excited states will be considered below in Section |V]) To this end, for each of the isospectral 
hamiltonians, we will calculate the neutron flux, assuming for simplicity that the states have 
zero widths and thus no interaction with the absorber is included. Due to these simplifying 
assumptions, a direct comparison with the existing experimental data is problematic, since 
an absorber is an important part of the actual experimental setup. Ideally, our results should 
be compared with a neutron flux measured without an absorber. Such an experiment would 
directly study properties of gravity itself, but a discussion of its feasibility lies beyond the 
scope of this paper. With this caveat, we will nevertheless proceed and use the data from 
Refs. treating it as an approximate guide to study our calculated fluxes. 



As explained in Section HTT| the ground state turning point z e * v G [10.3, 14.1] /im allows 
the potentials with p > 1.08 and p G [—1.13, —1.04] (see Fig. [p. Representing the three- 
dimensional wave function of a neutron as a product of a plane wave along the horizontal 
direction and the one-dimensional ground-state wave function ip\ (z; p) along the vertical (z) 
direction, the counting rate of the neutrons in the detector is proportional to the average 
horizontal velocity of the neutrons and to the integral 



We will call Di(z;p) the ground-state "density integral". More precisely, it is the probability 
of finding a neutron in the ground state ipi(z\p) above the mirror at any height less than 
or equal to z. An analogue of Di(z\p) is called "the neutron flux" in Ref. 0. Since, apart 
from normalization factors, all these quantities have the same z-dependence as the measured 
neutron count rate, we will use the terms "neutron flux", "count rate" and "density integral" 
interchangeably. 

We are interested in the behaviour of the density integral for the parameter p in the 
intervals compatible with the measured z\ v as described above, i. e. p > 1.08 and 
p G [—1.13,-1.04], and, for brevity, we will consider only the limiting values: p = 
1.08,-1.04,-1.13, as well as the flux calculated for the standard mgz potential (corre- 
sponding to p — > ±oo, as explained in Section [II])- The quantity plotted in Fig. [2] is a renor- 
malized density integral (z;p) = R\Di{z;p), where the constant R\ does not change the 



z 




(11) 
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FIG. 2: The ground-state density integral Eq. (Ilip . renormalized as explained in the text to 
facilitate the comparison with the measured neutron count rate [2]. The parameters p limit the 
range of the isospectral potentials compatible with experimental data for the ground state, as 
described in Fig. [TJ 



^-dependence of the flux and is used for simply to facilitate the comparison with the data 
from Ref. For completeness, we list the renormalization constants in the Appendix. 

Focusing on the behaviour of Di(z; p) near the lowest turning point z\ v G [10.3, 14.1] /zra, 
where the experimental flux exhibits a sharp step-like rise, we see that the rate of increase of 
the flux in the vicinity of z\ v strongly depends on the value of parameter p. This suggests 
that, within our model, we can select the most likely isospectral potential by comparing 
the corresponding density integral Di(z;p) with the shape of an experimental flux. A fairly 
good agreement with the calculated flux is obtained for the isospectral hamiltonians with 
— 1.13 < p < —1.04. It is noteworthy that, apart from the trivial data renormalization, 
at this stage we do not have any adjustable parameter as the range of the p was fixed 
independently in Fig. [TJ The description of the data in Fig. [2] is not perfect, due to the 
difficulties of a direct use of the existing data in our model, as explained in the beginning 
of this section. Note that D\(z\p = 1.08) cannot satisfactorily describe the data, even 
with a variety of different normalizations applied. The flux for the usual mgz potential 
also disagrees with the data; this disagreement emphasizes the importance of including an 
absorber in theoretical models based on the unmodified Newton potential, as was done in 
Refs. 0,0. 

We can see that the theoretical fluxes for —1.13 < p < —1.04 provide a good description 
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of the onset of the flux "step", which includes the data in the range of small z up to the 
point of the sharp increase at z ~ 13 fim. The agreement deteriorates for z > 16 fim, where 
effects of excited states, not included in Fig. [2j become important. The excited states are 
described by the wave functions Eq. (1101) . and in principle, many of them can be included 
in the flux. However, since only the first excited state has been observed [l|, 0] with an 
acceptable resolution, we will study its contribution in some detail. 

V. CONTRIBUTION OF THE EXCITED STATE 

In this section we consider the isospectral potentials allowed by both the ground and 
first excited state, i. e. those with p = 1.08 and p = —1.13, as explained in Section HTT1 
Assuming that a neutron is in a coherent superposition of the ground state ipi^p), given 
by Eq. ([8]), and the first excited state ipzfap), given by Eq. (ITU1) . its wave function can be 
written as 

ihafap) = ciipi(z;p) + c 2 ip 2 (z;p) , (12) 

where ci j2 coefficients describing the relative contributions of the states. By analogy 
with Eq. (TIT]) , we consider the "two-state density integral" 

2 

D 12 (z;p) = Jdz'tf 2 (z'-p), (13) 
o 

which implicitly depends on ci )2 . 

We will explore to what extent the two-state density integral Di 2 (z;p) may be used to 
describe of the experimental flux up to the second turning point z^ v ■ Since the neutron 
populations in separate states were not directly measured [2] , one can treat C\ and c 2 as free 
parameters. However, taking into account the simplicity of our model, we will not attempt 
to find the best quantitative fit to the data. Rather, our aim is to illustrate the contribution 
of the excited state for several contrasting choices of C\ and c 2 . Similar to the ground-state 
density integral Di(z;p), we renormalize Di 2 (z;p) for comparison with the data and thus 
introduce Df 2 (z;p) = Ri 2 D\ 2 (z; p). As explained in Section HV] this renormalization does 
not affect the important ^-dependence of the flux. The method of fixing the renormalization 
constants R\ 2 and their values are described in the Appendix. 

Fig. [3] shows the renormalized two-state density integral D± 2 (z;p) calculated for the 
isospectral potentials with p = 1.08 and p = —1.13, as well as for the mgz potential, 
and for several sets of Ci j2 . The bottom panels correspond to the cases where the excited 
state contributes with the same weight as the ground state, with two selections of signs: 
ci = ±c 2 . The excited state contribution is emphasized in the middle panels (|c 2 /ci| = 2), 
and attenuated in the bottom panels (|c 2 /ci| = 1/2). When c\ and c 2 have opposite signs, 
the density integrals of the isospectral potential with p = —1.13 lie closer to the data than 
those of the potential with p = 1.08 or the Newton potential mgz. For c\ and c 2 of the 
same sign, the potential with p = 1.08 or the mgz perform better than that with p = —1.13. 
Among all calculations shown in Fig. [3], the case with p = —1.13 and C\ = — c 2 = 1 seems to 
describe the data better than the others. However, the agreement of the calculated fluxes 
with the data over the whole range of heights up to z ~ 30 fim is rather poor. Some of the 
reasons for this discrepancy have been already discussed in Section HVl such as the absence 
of an absorber in the present model. On the other hand, an attractive feature of the model is 
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FIG. 3: The renormalized two-state density integral Eq. (|13p for the isospectral potentials with 
p = 1.08 (left panels) and p = —1.13 (right panels), compatible with the measured turning points 
as shown in Fig. [TJ The fluxes corresponding to several sets of constants ci^ are shown 
from the top to bottom panels, to illustrate several distinct superpositions of the ground and first 
excited states. The data points are from Ref. 0], and the flux renormalization factors are listed in 
Appendix. 



that it uses only two parameters {c\p) to incorporate both the ground and the first excited 
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state. Such a simplified approach has the advantage of dealing exclusively with neutron 
properties under gravity, circumventing a possible explanation of the experimental results 
of [l| as a "neutron in a box" quantum system without any potential (see discussion in 
Refs. 0, El])- Assuming that the neutrons are subjected to gravity alone, without the influ- 
ence of an absorber, Figs. and ([3]) show that the fluxes calculated with the isospectral 
potentials can provide a qualitative description of the essential features of the data. 

Contributions of the higher excited states can be incorporated analogously by calculating 
a multi-state density integral Di 2 ...n (for N states), with iV adjustable parameters c\ y 2,...,N- 
However, since the experimental resolution of the higher excited states is at present quite 
poor [2], we will not extend the present analysis beyond the first excited state. 



VI. DISCUSSION AND OUTLOOK 

There may be two possible interpretations of the results presented in this paper. One 
possibility is to regarded this model as an effective approach which mimics models such as 
those of Refs. 0, 0], in which neutrons move under the usual Newton force, where a bottom 
mirror and a top absorber would be described in detail. Another possibility is to view the 
isospectral potentials considered here as true modifications of the Newton gravitational law 
at a micrometer scale. The characteristic length over which the deviation from the linear 
potential is observed is given by / = [h 2 /(2m 2 g)] lj ' 3 , which scales as ~ mr 2 / 3 . This means 
that in experiments involving objects heavier than neutrons, modifications from the Newton 
law would be present on a shorter scale, which could explain why such modifications are 
difficult to observe for atoms or macro-objects. 

It is interesting to note that in the vicinity of the turning point z\ ~ 13/xm the isospectral 
potentials with —1.13 < p < —1.04 are alternately more attractive or more repulsive than 
the unmodified Newton potential mgz. It is argued 0, H[ that both the geometry of extra 
dimensions and the nature of quantum fields may influence whether the modification is 
attractive or repulsive. 

Our results suggest that the isospectral potential which most closely describe the prop- 
erties of the ground the first excited states, including the turning heights and (only qual- 
itatively)the flux profile, is that with the parameter p pa —1.13. A peculiar feature of 
this potential is a barrier of a certain width between z = and the lowest turning point 
z\ pa 13/im (see the right panel of Fig. [I]), which means that the small probability of finding 
a particle near z ~ is a purely quantum effect of barrier penetration. Although such a 
gravitational potential looks rather unusual from the classical mechanics point of view, it is 
perfectly legitimate since it has the same energy levels as the Newton potential and is also 
consistent with the experimentally observed properties of the neutron gravitational ground 
state. As mentioned above, it would be interesting to carry out an experiment similar to 
that of Refs. [l|, 0] , but measuring the vertical distribution of the neutrons in a setup without 
an absorber. This would provide a more direct test of the present model. 

One recently proposed [19J direction of studying the excited gravitational states is to 
measure fluxes of spin-polarized ultra-cold neutrons in the presence of magnetic and gravi- 
tational fields, using the high-precision GRANIT spectrometer j20[. Another improvement 
can be due to the use of new sources of ultra-cold neutrons, notably the currently developed 
UCN source at TRIUMF [2lj], based on the earlier spallation design described in Ref. (22|. 
This and other high-density sources are expected to yield data with smaller error bars and 
a better resolution of the states. 
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APPENDIX A: RENORMALIZATION CONSTANTS FOR THE DENSITY IN- 
TEGRALS 



TABLE I: Renormalization constants for comparing the density integrals with data. 



Renormalization factor 


"Best fit" flux 


#1 = 0.0043 


p = -1.13 


R 12 = 0.0071 (ci = 1,Q2 = 1) 


p = ±oo(mgz) 


R 12 = 0.0048 (ci = 1,Q2 = 2) 


p = 1.08 


R12 = 0.0125 (ci = l,c 2 = 1/2) 


p = ±oo(mgz) 


#12 = 0.0096 (c! = l,c 2 = -1) 


p= -1.13 


#12 = 0.0042 (ci = l,c 2 = -2) 


p= -1.13 


#12 = 0.0135 (ci = l,c 2 = -1/2) 


p= -1.13 



In Table [I] we summarize the constants R\ and R12 which are used to renormalize the 
density integrals Eqs. (fill and (f!3l) for comparison with the data in Figs. (ED and (ED, respec- 
tively. In practice, the renormalization constants were chosen so that, for each combination 
of ci )2 , one of the calculated fluxes described the data in the best achievable way; these 
"best fits" are listed in the right column. After that, the same values of R\ and -R12 were 
used for all values of the parameter p. It should be stressed that we did not attempt to do 
a quantitative fit to the data over the whole range of allowed p and ci^, which is justified 
within the present simplified model. 
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